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Abstract 

The  buckling  loads  of  8-ply  unstiffened  laminated  (graphite-epoxy) 
cylindrical  panels  were  determined  experimentally.  The  results  were 
compared  to  the  linear  bifurcation  loads  as  determined  by  the  STA6SC-1 
Computer  Code.  The  analysis  Included  three  different  ply  orientations, 
five  different  aspect  ratios  (length^chord),  and  one  set  of  boundary 
conditions;  clamped  along  the  top  and  bottom  curved  edges  and  simply 
supported  along  the  vertical  sides.  The  experimental  tests  were  con¬ 
ducted  by  the  Air  Force  Flight  Dynamics  Laboratory  using  a  redesigned 
test  set-up.  Fairly  good  agreement  was  obtained  between  the  analytical 
and  experimental  loads,  with  the  linear  bifurcation  values  being  3%  to 
32%  higher  than  the  experimental.  It  was  evident  that  the  buckling 
pattern  observed  experimentally  did  not  In  all  cases  correspond  to  the 
symnetrical  pattern  generated  by  STAGSC-1.  This  was  due  mainly  to  the 
test  fixtures  Inability  to  evenly  distribute  the  stress  along  the 
simple  support  mechanism. 
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ANALYTICAL/ EXPERIMENTAL  LINEAR  BIFURCATION 
OF  CURVED  CYLINDRICAL  COMPOSITE  PANELS 

I.  INTRODUCTION 

Background 

Composite  materials'  high  strength  to  weight  ratio  make  it  par¬ 
ticularly  well  suited  for  aircraft  structural  use.  The  USAF's  forward 
swept  aircraft  design  [1]  is  an  excellent  example  of  the  materials' 
advantage  in  that  through  proper  composite  "tailoring",  the  problem  of 
aeroelastic  divergence  can  be  practically  avoided  without  the  enormous 
weight  gain  associated  with  conventional  metals.  However,  in  comparison 
with  metals  the  science  of  composites  is  still  in  its  infancy,  therefore, 
engineers  are  reluctant  to  base  their  design  solely  on  computer  generated 
data  without  some  comparison  to  realistic  emperical  data. 

Numerous  articles  have  dealt  with  the  analysis  of  composite  cylind¬ 
ers  and  flat  plates  [2-4].  But  due  to  manufacturing  difficulties  and 
the  lack  of  experimental  test  fixtures,  literature  on  actual  experimenta¬ 
tion  of  cylindrical  composite  panels  is  relatively  scarce.  One  of  the 
first  and  most  influential  articles  on  the  subject  was  presented  by 
Wilkens  [5]  in  1974.  Using  a  unique  test  set-up,  he  determined  the 
buckling  loads  of  different  ply  orientations  associated  with  one  width 
to  length  ratio  (referred  to  as  the  aspect  ratio).  His  results  show 

percent  knockdow,  factors  ((1  -  thgoI^^Ml’bSckllnS^liS’^)  *  ”"9- 

ing  from  1%  to  66%.  He  also  found  that  multidirectional  laminates 
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retain  a  high  percentage  of  their  buckling  strength  after  repeated 
buckling.  More  recently,  Bauld  [6]  used  a  modification  of  Wil kens' 
test  fixture  to  study  the  experimental  buckling  loads  of  composite 
cylindrical  panels  having  two  different  ply  orientations  and  five 
different  aspect  ratios.  In  particular,  he  found  that  circumferential 
imperfections  in  a  panel  serve  to  lower  its  resistance  to  buckling. 

This  thesis  is  a  continuation  of  the  work  done  by  Becker  [7],  with  a 
redesigned  experimental  test  fixture  which  allows  an  even  distribution 
of  applied  load  and  thus  a  more  reliable  set  of  experimental  results. 
The  work  reported  herein  involved  a  new  set  of  aspect  ratios  allowing 
for  a  broader  appreciation  of  panel  buckling. 

Purpose 

The  purpose  of  this  thesis  is  to  compare  experimentally  determined 
buckling  loads  of  cylindrical  composite  panels  to  the  linear  bifurca¬ 
tion  loads  calculated  by  the  STAGSC-1  Computer  Code.  The  effect  of 
ply  orientation  and  aspect  ratio  will  be  the  two  major  variables.  The 
results  will  be  used  to  further  the  understanding  of  composite  buckling 
behavior  when  subjected  to  design  loads  [8]. 


Scope 

A  total  of  30  panels  were  tested  and  aspect  ratios  of  ^  (8"xl6"), 

^  (12"xl6"),  1  (16"xl6"),  ^  (16"xl2"),  and  2  (16"x8")  were  physically 
included.  All  panels  consisted  of  8-ply  graphite-epoxy  with  experiment¬ 
ally  determined  material  properties  of  [7]: 


Ej  =  20.5  X  10° 

Eg  =  1.3  X  10® 

G^2  =  0.75  X  10® 

Vi2  *  0.335 

The  three  different  ply  orientations  used  are  (±45)25,  (90,0)25,  and 
(90,±45,0)5.  Figure  1  shows  a  typical  panel  and  defines  the  panel 
notations  and  sign  conventions.  The  panel  boundary  conditions  incor¬ 
porated  in  this  work  consists  of  a  clamped  top  curved  edge  (u=free, 
v=w=w  =0),  a  clamped  bottom  curved  edge  (u=v=w=w  ^=0),  and  simple 

supports  along  the  vertical  sides  (u=v=w  =free,w=0).  See  Fig.  1 

»y 

for  positive  displacement  directions.  It  should  be  noted  that  the  30 
panels  represent  a  duplication  of  each  unique  ply  orientation  and  aspect 
ratio,  thus  an  average  experimental  buckling  load  can  be  obtained  for  a 
given  set  of  parameters. 


t  =  thickness  =  8  plies  @  0.005"  =  0.04" 
R  »  radius  =  12" 

C  =  chord  length 
L  ■  length 

x,y,z  =  structural  coordinate  directions 
9  s  ply  orientation 
u,v,w  =  displacements 
L/C  =  aspect  ratio 


Figure  1 


Panel  Notation 


II.  THEORY 


Bifurcation  Buckling 

Bifurcation  involves  the  transfer  of  membrane  strain  energy  to 
bending  strain  energy.  Bushnell  [9]  has  shown  that  a  thin  shell 
(cylindrical  panel  in  this  thesis)  can  absorb  a  great  deal  of  membrane 
strain  energy  with  little  deformation.  The  shell  must  deform  much  more 
for  it  to  absorb  an  equivalent  amount  of  bending  strain  energy.  It  is 
through  these  larger  deformations  that  the  transfer  of  energies  will 
occur  resulting  in  bifurcation. 

The  above  explanation  is  best  illustrated  by  relating  it  to  a 
panel's  load-displacement  curve,  or  synonymously,  its  equilibrium  path 
[10].  Figure  2  shows  the  load  versus  end  displacement  relationship  for 
a  thin-walled  cylindrical  panel  subjected  to  axial  compression.  The 
primary  equilibrium  path  is  intersected  by  a  secondary  path  which  is 
characterized  by  a  dramatic  drop-off  in  load  and  a  negative  slope.  The 
intersection  is  the  bifurcation  point  and  represents  the  buckled  config¬ 
uration  of  the  panel.  At  bifurcation,  the  shell's  equilibrium  equations 
have  multiple  solutions.  The  so-called  secondary  path  represents  a 
stable  solution  while  the  portion  of  the  load-displacement  curve  that 
lies  above  the  bifurcation  point,  referred  to  as  the  primary  path,  is 
unstable. 


x,y,z  =  Structural  Axes 
1,2,3  =  Lamina  Principal  Axes 


Figure  3 


Definition  of  Coordinate  Systems  [7] 


In  any  other  coordinate  system  with  the  fiber  axis  oriented  at 
some  angle,  e,  with  respect  to  the  structural  axes,  the  stresses  are 


where  are  the  transformed  reduced  stiffness 

Q  =  Q  cos'*0  +  2(Q  +  2Q  )  sln^e  cos^e  +  Q  s1n'*0 

11  11  12  66  22 

0  =  (Q  +  Q  -  4Q  )  s1n*0  cos^0  +  Q  (s1n'*0  +  cos'*0) 

12  11  22  66  12 

Q  *  Q  s1n'*0  +  2(Q  +  2Q  )  s1n*0  cos*0  +  Q  cos'*0 

22  11  12  66  22 

Q  =  (Q  -  Q  -  2q  )  Sln0  cos*0  +  (Q  -  Q  +  2Q  ) 

16  11  12  66  12  22  66 

s1n^0  COS0 

Q  =  (Q  -  Q  -  2q  )  sin®0  cos0  +  (q  -  q  +  2q  ) 

26  11  12  66  12  22  66 

s1n0  cos®0 

Q  *  (q  +  q  -  2q  -  2q  )  sin^0  cos*0  +  q  (sin'*0 


(! 


+  COS**0) 


Equation  (4)  can  be  thought  of  as  the  stress-strain  relations  for 
the  K**'  Layer  of  a  multilayered  laminate.  Thus,  it  can  be  written  as 

tol^  =  CQ]^  {e},^  (6) 

To  expand  the  discussion  to  a  multilayered  laminate,  a  few  assump¬ 
tions  must  be  made.  First,  the  laminate  is  presumed  to  consist  of 
perfectly  bonded  laminae  so  that  no  lamina  can  slip  relative  to  another. 
Secondly,  adherence  to  the  Kirchhoff-Love  hypothesis  dictates  that 
sections  normal  to  the  panel's  mid-surface  remain  plane  and  normal  to 
that  surface  after  bending  [11].  These  assumptions  lead  to  the  follow¬ 
ing  strain-curvature  relationship  for  a  laminate 


where  the  o  superscript  indicates  the  mid-surface  strains,  the  k's  are 
mid-surface  curvatures  and  Z  represents  a  distance  from  the  panel  mid¬ 
surface.  For  a  cylindrical  panel  with  moderately  large  displacements 
and  also  moderately  large  rotations  of  tangents  to  its  mid-surface,  the 
above  strains  and  curvatures  are  given  by  Sanders'  kinematic  relations 
without  initial  imperfections  [13] 


U  +  4  +  i  4^ 

,x  .  2  ’x  2 


,X 


+  u 


.y 


+  4  4 
^x 


(8 


+ 


1 

r 


where  u,v,w  denote  the  axial,  circumferential,  and  radial  components 
of  displacement,  respectively,  of  the  panel's  mid-surface.  The  4's 
are  rotation  components  and  are  expressed  in  terms  of  the  displacements 
as 


♦y  ■  -".y  *  7  (9) 

*  *  7  <'',x'“,y> 


where  r  is  the  panel's  radius  of  curvature. 

Substituting  Eq.  (7)  into  Eq.  (4),  the  stresses  in  the  Layer 
can  be  expressed  in  terms  of  the  laminate  middle  surface  strains  and 
curvatures  as 


m 

X 

>  * 

y 

^^xy 

s 

c  L 

11  12  16 


12  22  26 


16  26  66 


Uk  I 


fe"  n 

r  1C  ’  'I 

®X 

’'x 

e® 

Z  i 

)C 

1 

y 

y 

Y® 

''^xy^ 

.2«xy. 

(10) 
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Now  that  the  stresses  in  each  layer  are  known,  the  resultant  forces 
and  moments  acting  on  a  laminate  are  obtained  by  integration  of  the 
stresses  in  each  lamina  through  the  laminate  thickness.  For  example. 


and 


ftl2 

»  a  ZdZ  (12) 

X  J_t/2  ^ 


and  is  a  force  and  moment  per  unit  length  (width)  of  the  cross- 
section  of  the  laminate  as  shown  in  Fig.  4,  and  are  defined  as 


where  N  is  the  total  number  of  lamina,  and  Z.  and  Z.  ,  are  defined  in 


Figure  4 


Forces  and  Moments  on  a  Laminate 


Figure  5  Geometry  of  an  N-Layered  Laminate 


The  integration  indicated  by  Eq's.  (13)  and  (14)  can  be  rearranged 

to  take  advantage  of  the  fact  that  the  stiffnesses  matrix,  (J. for 

*  0 

each  layer  is  constant  within  that  layer.  Thus,  the  stiffness  matrix 
goes  outside  the  integration  over  each  layer  but  is  within  the  summa¬ 
tion  of  force  and  moment  resultants  for  each  layer.  Also,  recall  that 

Ey*  e®,  y“  »  Ky.  and  K  are  not  functions  of  Z  but  are  middle 
X  y  xy  X  y  xy 

surface  values  so  they  too  can  be  removed  from  under  the  summation  signs. 
Thus,  the  force  and  moment  equations  become 


14 


where 


N 


The  A.,  are  called  extensional  stiffnesses,  the  B..  are  called  coupling 

I  J  I  J 

stiffnesses,  and  the  D. .  are  called  bending  stiffness.  The  mid-surface 

*  J 

strains  and  curvatures  are  expressed  in  terms  of  displacements  (Eq's. 
(8)  and  (9))  to  yield 


STAGSC-1  Theory 

STAGSC-1  (Structural  Analysis  of  General  Shells)  is  a  computer 
code  developed  by  the  Lockheed  Palo  Alto  Research  Laboratory  to  analyze 
general  shells  under  various  static,  thermal,  and  mechanical  loading 
[14].  The  approach  used  by  STAGSC-1  is  an  energy-based  finite  element 
analysis.  According  to  the  energy  method,  a  system's  total  potential 
energy  is  used  to  derive  its  equilibrium  equations  from  which  stability 
can  be  determined  by  the  solution  to  an  eigenvalue  problem. 
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A  shell  s  total  potential  energy  is  equal  to  its  internal  strain 
energy  minus  the  product  of  the  external  forces  and  their  respective 
deflections. 

An  element's  strain  energy  is  given  by  [13]: 

^  Area  °  ° 

where 


The  expressions  for  the  mid-surface  strains  are  given  in  Eq's.  (19)  and 
[N]  was  presented  in  Eq.  (17).  In  general,  the  strain  vector  [e]„  is 


a  function  of  the  mid-surface  displacements  u,v,w,  the  first  order  de¬ 
rivatives  of  u,v,w  and  the  second  order  derivatives  of  w.  Therefore, 
we  can  let  a  vector  [d]  represent  this  functional  dependence  of  [e]^ 
on  the  displacements  by 


By  using  vector  [d],  Bauld  [13]  carried  out  the  integration  of  Eq.  (20) 
and  found  that  the  expression  for  strain  energy  is  comprised  of  three 
distinct  parts.  The  first  part  is  quadratic  in  displacements,  the  sec¬ 
ond  part  is  cubic  in  displacements,  and  the  third  part  is  quartic  in 
displacements.  In  a  finite  element  analysis,  vector  [d]  is  represented 
by  an  element's  shape  functions  and  nodal  degrees  of  freedom.  Thus,  the 
strain  energy  can  be  expressed  in  terms  of  an  element's  degrees  of  free¬ 
dom.  An  element's  external  potential  energy  is  obtained  by  using  a 
similar  finite  element  analysis  on  its  external  forces.  The  strain 
energy  and  external  potential  energy  is  combined  to  give  the  total 
system  potential  energy,  V,  and  can  be  written  in  the  form  given  by 
Bauld  [13]  as 

*  '  ^  '  %'>s  <22) 

A  is  the  system's  linear  stiffness  matrix  with  no  dependence  on  the 

I  d 

displacement  vector,  q.  Nl^^  and  N2^g  are  matrices  with  linear  and 
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quadratic  dependence,  respectively,  on  displacement.  Is  the  surface 
force  vector. 

The  principle  of  total  potential  energy  states  that  the  equilibrium 
configuration  of  a  conservative  mechanical  system  corresponds  to  a 
stationary  value  of  the  total  potential  energy  of  the  system  [15]. 
Therefore,  taking  the  first  variation  of  Eq.  (22)  and  setting  It  to 
zero,  we  have  a  set  of  nonlinear,  algebraic  equations  of  the  form 

(A„  ^  .  0  (23) 

Loss  of  stability  (bifurcation),  results  when  the  second  variation  of 
the  system's  total  potential  energy  ceases  to  be  positive  definite,  or 

DET  (A^5  +  =  0  (24) 

Equation  (24)  Is  used  by  STAGSC-l  to  solve  the  eigenvalue  problem  of 
the  form  [10] 

[A]  +  x[B]  +  X^C]  =  0  (25) 

where  [B]  and  [C]  represent  nonlinear  stiffness  matrices  In  unknown 
displacements  and  products  of  displacements,  respectively.  For  a 
linear  analysis,  x.  Is  the  proportionality  constant  of  a  convenient 
load  level  used  In  the  equilibrium  Eq's.  (23)  to  solve  for  the  unknown 
displacements.  The  quantities  A,  B,  C,  or  equivalently,  A^^,  Nl^^,  and 
N2^g  are  calculated  once  based  on  the  equilibrium  displacements.  Fin¬ 
ally,  the  load  proportionality  parameter,  x.  Is  Incremented  until  a 
sign  change  on  the  left  side  of  Eq.  (24)  occurs,  signifying  bifurcation. 


1 


p 
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III.  FINITE  ELEMENT  MODELING 

Boundary  Conditions 

The  STAGSC-1  computer  code  has  a  large  variety  of  panel  boundary 
conditions  to  choose  from  ranging  from  the  most  rigid  to  the  most 
flexible  conditions  that  can  exist  in  an  experimental  test.  Reference 
[7]  gives  analytical  results  associated  with  the  full  range.  The 
author's  choice  of  a  simple  support  along  the  vertical  sides  and 
clamped  along  the  curved  edges  was  strictly  made  to  approximate  the 
conditions  used  experimentally. 

Grid  Size 

The  determination  of  a  finite  element  grid  size  is  crucial  to 
obtaining  accurate  results.  Nelson  [16]  found  that  to  obtain  accurate 
results  there  should  be  at  least  five  node  points  per  each  half  sine 
wave  of  the  buckled  pattern  in  the  circumferential  direction.  Becker 
[7]  used  this  information  and  found  that  a  0.5  x  0.5  inch  grid  repre¬ 
sented  an  ideal  trade-off  between  accuracy  and  computer  economics.  The 
0.5  X  0.5  inch  grid  is  used  in  this  thesis,  see  Fig.  6.  As  a  means  of 
comparison,  the  bifurcation  load  of  the  (±45)2^,  ^  aspect  ratio  panel 
was  analytically  determined  by  using  a  1.0  x  1.0  inch  grid  and  a  0.5  x 
0.5  inch  grid.  The  bifurcation  loads  differed  by  15%.  However,  the 
finer  grid  required  over  six  times  the  computer  time  at  10  times  the 
cost.  See  Appendix  B  for  a  sample  STAGSC-1  computer  program. 
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Finite  Element  Selection 


Figure  7  shows  the  410  Quadrilateral  Plate  (QUAF)  element  used  In 
this  thesis.  The  QUAF  element  was  specifically  developed  to  eliminate 
the  Interelement  displacement  Incompatibility  that  exists  when  flat 
elements  are  used  for  curved  surfaces  [17].  The  410  element  accomplishes 
this  by  Increasing  the  order  of  the  polynomial  expression  for  the  Inplane 
Displacement  Field  by  virtue  of  a  Normal  Rotation  Component,  e^.  The 
result  Is  that  each  element  has  6  degrees  of  freedom  (3  translation  and 
3  rotation)  per  node  for  a  total  of  24.  The  410  element  works  well  for 
a  linear  analysis  but  does  present  problems  when  nonlinearity  Is  being 
studied. 

For  nonlinear  problems,  the  QUAF  411  element.  Fig.  8,  Is  the  better 
choice  because  of  the  y  degrees  of  freedom,  which  represent  the  normal 
rotation  of  the  two  Intersecting  grid  lines.  A  comparison  of  linear  bi¬ 
furcation  loads  was  made  using  the  two  elements.  It  was  found  that  the 
load  varied  4%  for  the  (±45)25,  2-aspect  ratio  panel,  and  only  .23%  for 
the  (90,0)25,  2-aspect  ratio  panel.  However,  the  411  element  was  Zj 
times  more  expensive.  The  small  difference  In  loads  correspond  to  the 
results  of  Snead  [18]  who  showed  that  the  difference  In  loads  associat¬ 
ed  with  the  two  elements  Is  dependent  upon  panel  ply  orientation  and 
the  number  of  elements  In  the  finite  element  grid.  In  general,  the 
greater  the  number  of  elements,  the  smaller  the  difference. 


Figure  7 
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IV.  EXPERIMENTAL  PROCEDURE 


Test  Panels 

All  of  the  panels  tested  were  made  from  Narmco  T300/5208  Graphite 
Epoxy  and  manufactured  in  1979.  They  were  originally  laid  up  into  17" 

X  36"  panels  which  were  cut  to  the  desired  size  by  means  of  a  specific¬ 
ally  jigged  radial  arm  saw  with  a  diamond  tip  blade.  No  special  tech¬ 
niques  were  used  to  insure  "perfect"  panels  for  the  test,  although  every 
effort  was  made  to  fabricate  and  cut  the  panels  as  close  as  possible  to 
the  desired  test  configuration.  Therefore,  the  normal  dimensional  var¬ 
iations  (thickness,  length,  and  width)  inherent  in  all  production  proc¬ 
esses  were  present  [7].  Another  variable  that  probably  had  some  effect 
upon  the  results  was  the  age  of  the  panels.  Although  they  were  stored 
in  controlled  humidity  and  temperature  conditions,  these  factors  should 
be  considered  in  critical  design  applications.  Environmental  variables 
are  beyond  the  scope  of  this  thesis. 

Test  Set-Up 

In  order  to  vary  the  boundary  conditions  and  aspect  ratios,  an 
experimental  device  needed  to  be  designed.  The  test  fixture  developed 
is  a  minor  modification  of  the  General  Dynamics  original  1974  design  [5] 
The  device  is  shown  on  Fig.  9.  The  reader  can  observe  the  capability  of 
this  fixture  to  allow  clamping  along  the  panel's  curved  edges  by  virtue 
of  the  small  blocks  shown  in  Fig.  9.  These  blocks  were  forced  against 
the  panel  by  a  series  of  set  screws.  The  set  screws  were  sufficiently 


tightened  by  a  screwdriver  to  insure  no  panel  movement.  The  metal 
flange  located  at  the  center  of  the  base  plate.  Fig.  9,  was  used  to 
insure  exact  centering  of  the  plate  with  respect  to  a  duplicate  top 
plate  (for  clamping  along  the  panel's  top  curved  edge).  This  technique 
avoids  load  eccentricities.  A  simple  support  condition  along  the  panel's 
vertical  sides  is  obtained  by  the  knife-like  edges  shown  in  Fig.  9  and 
10.  A  panel  was  inserted  between  the  two  knife  edges  and  then  the  set 
screws  were  finger  tightened  until  the  knife  edges  were  just  touching 
the  panel.  This  device  prevented  panel  movement  in  the  w  direction  but 
allowed  movement  in  the  u  and  v  directions.  In  all  tests  the  panels 
were  clamped  to  the  bottom  and  top  plates  first,  then  the  simple  supports 
were  installed  and  set  screws  finger  tightened.  Figure  11  shows  the 
complete  set-up  for  the  experimental  buckling  of  a  curved  cylindrical 
panel.  The  panel  support  system,  with  panel  installed,  is  located  in 
the  center  portion  of  the  set-up.  The  load  is  applied  through  a  hydraul¬ 
ic  ram  to  the  lower  steel  bed  which  moves  up  the  vertical  sleeves  of  the 
machine.  The  i  inch  panel  protrusion  at  the  top  of  the  simple  supports, 
see  Fig.  12,  allows  the  panel  to  be  compressed  against  the  top  hori¬ 
zontal  cross-head.  The  observed  panel  buckling  pattern  associated  with 
the  existing  set-up  was  not  symmetrical  as  was  previously  assumed.  This 
was  due  to  the  fact  that  the  i  inch  panel  free  space  was  at  the  top  of 
the  simple  supports.  It  was  found  that  when  the  i  inch  space  was  at  the 
bottom  of  the  support  mechanism  (obtained  by  repositioning  the  knife 
edges),  the  desired  buckling  symmetry  did  occur. 


One  might  observe  the  vertical  LVDT  (Linear  Variable  Differential 
Transducer)  located  on  the  right  side  of  the  set-up  in  Fig.  11,  used 
to  measure  panel  end-shortening.  Furthermore,  a  matrix  of  LVDT's  have 
been  placed  horizontally  on  the  back  of  the  panel  as  shown  in  Fig.  13, 
to  measure  radial  displacements.  Additionally,  6  axial  strain  gages 
(back-to-back  in  3  locations)  were  attached  to  each  of  the  30  panels. 
Included  within  the  experimentation  is: 

1.  determination  of  buckling  load, 

2.  strain  measurements  on  inner  and  outer  faces  of 
the  panel  at  discretely  located  positions, 

3.  radial  displacement  measurements. 

These  three  measurements  were  incorporated  into  three  techniques  for 
determining  the  buckling  load 

1.  load  vs  end-shortening, 

2.  strain  reversal, 

3.  Southwell  Plot  (see  Appendix  A). 

Each  set  of  measurements  were  carried  out  for  all  30  panels.  The 
data  was  stored  in  a  CEL  86  Computer  in  a  real  time  mode  and  printed 
out  later.  The  data  was  also  displayed  on  a  computer  terminal  CRT  so 
that  the  load  and  displacements  could  be  viewed  during  the  testing. 

Test  Procedure 

Once  the  panel  was  installed,  the  instrumentation  (LVDT's  and 


strain  gages)  was  zeroed  out  electronically  and  the  test  begun.  The 
load  was  applied  and  the  data  was  sampled  two  times  per  second  for 
most  of  the  panels.  When  the  first  major  buckle  occurred,  the  location 
was  noted  and  the  buckled  pattern  photographed.  The  loading  was  con¬ 
tinued  In  order  to  get  some  portion  of  the  postbuckling  curve.  When 
the  top  plate  contacted  the  side  supports,  or  when  panel  breakage  made 
continued  loading  Impractical,  the  load  was  removed. 


V.  RESULTS 


Panel  Identification 

The  panels  used  in  this  thesis  are  identified  in  the  following 
manner: 

1  1  or  Duplicate  Panel  2 
Orientation 
=  (90.±45.0)g 

=  (90.0)25 
=  (±45)25 

1  Size  (length  x  chord) 
19=16"xl6";  AR»1 
29=16"x8'';  AR«2 
39=8''xl6";  AR=l/2 
49=16''xl2”;  AR*4/3 
59=12"xl6";  AR=3/4 


A  three  digit  number  is  used  when  the  average  results  of  both 
duplicate  panels  are  referenced.  As  an  example.  392  means  the  average 
buckling  load  or  end-shortening  deflection  of  panels  3921  and  3922  are 
being  referenced. 


Exper i menta 1 /Ana 1 y t i ca 1  Compar 1  son 


Table  I  shows  the  experimental  and  analytical  results  for  all  30 


panels  tested.  Table  II  Is  a  summary  of  Table  I,  showing  only  the 
average  knockdown  factors  associated  with  a11  three  methods  used  to 
determine  the  buckling  loads.  The  percent  knockdown  factor  used  In  this 
thesis  is  defies.  ,s  (1  -  *  ‘“O-  end- 

shortening  technique  (E-S).  was  the  most  reliable  source  of  experimental 
data  with  knockdown  factors  ranging  from  2.9%  to  31.8%.  The  buckling 
load  was  determined  to  be  that  maximum  load  which  caused  the  buckling 
configuration  of  the  panel  and  a  distinct  drop-off  In  load  on  a  plot  of 
Its  load  vs  end-shortening  curve.  The  load  vs  end-shortening  curve  for 
each  panel  Is  shown  In  Fig's.  14  to  27.  From  these  curves  It  can  be 
noticed  that  the  panels  do  not  have  post-buckling  strength,  which  was 
similarly  reported  by  Becker  [7].  Panel  1911  was  not  Included  In  the 
analysis  because  panel  breakage  produced  such  scattered  load  vs  end¬ 
shortening  data  that  a  linearization  of  the  curve  proved  Impractical. 
Panels  2921,  2922,  3922,  and  5921  were  not  Included  In  the  calculation 
of  E-S  knockdown  factors  because  their  experimental  buckling  loads  were 
higher  than  the  analytical  values.  This  was  most  likely  due  to  panel 
cutting  Irregularities.  These  panels  were  cut  a  little  wider  than  the 
rest  of  the  panels  praducing  a  smaller  clearance  between  the  vertical 
edges  of  the  panels  and  the  back  of  the  simple  support  device  shown  In 
Fig.  10.  As  the  load  was  applied  the  panel's  freedom  of  motion  In  the 
Y-dIrectlon  was  surpressed  resulting  In  a  stiffen  panel.  A  new  set  of 
boundary  conditions  (with  v»0)  were  programmed  analytically.  The  re- 
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Table  I  Key 


Panel  Breakage  Procluded  Meaningful  Data 

Linear  Relationship  of  Ld  vs  Circumferential  Displ. 

Unreducable  Data 

Buckling  Load  Higher  Than  Analytical  Bifurcation  Load 
Strain  Gages  Not  Placed  Back-to-Back 
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Table  II  Key 


I 


,  ,  _  End-Shortening  Buckling  Load 

*  ^  "  Analytical  Bifurcation  Load 

2  1  -  Southwell  Buckling  Load 

*  ^  "  Analytical  Bifurcation  Load 

2  ,  Strain  Reversal  Buckling  Load 

*  ^  "  Analytical  Bifurcation  Load 


4. 

5. 

6. 

7. 

8. 
9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 


Panel  1911  not  included  due  to  breakage 

Both  2921  and  2922  had  higher  experimental  buckling  loads  than  analytical 
3922  not  included;  higher  than  analytical  load 
1922  not  included;  linear  data 
Linear  data  for  both  panels 


2921  not  included 

2922  not  included 
3932  not  included 
4922  not  included 
4922  not  included 
5922  not  included 
5921  not  included 
5921  not  included 


higher  than  analytical  load 
unreducable  data 
linear  data 

Southwell  load  higher  than  analytical 
unreducable  data 
linear  data 

strain  gages  not  placed  back-to-back 
higher  than  analytical  load 
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Figure  16  Analytical  and  Experimental  Results  for  193  Panel 


Analytical  and  Experimental  Results  for  291  Panel 


Analytical  and  Experimental  Results  for  3921  Panel 
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Figure  21  Analytical  and  Experimental  Results  for  393  Panel 
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Analytical  and  Experimental  Results  for  593  Panel 


suiting  analytical  bifurcation  loads  were,  in  fact,  higher  than  the 
experimental  buckling  loads.  This  indicates  that  the  actual  experiment¬ 
al  boundary  conditions  for  these  four  panels  did  not  correspond  to  the 
desired  simple  support  conditions.  The  knockdown  factors  associated 
with  the  new  boundary  conditions  (v=0)  were  11%  and  13.7%  for  the  292 
and  5921  panels,  respectively. 

The  strain  reversal  technique  (S-R)  provided  a  good  check  of  the 
E-S  method,  with  knockdown  factors  ranging  from  5.1%  to  41%.  However, 
S-R  was  particularly  susceptible  to  computer  channel  "noise"  which  on 
occasion  produced  unreducable  data.  Figure  28  shows  a  diagram  of  the 
strain  reversal  process.  As  the  load  is  increased,  the  back-to-back 
strain  gages  register  essentially  the  same  increase  in  strain  values. 
When  buckling  occurs,  the  strains  of  the  respective  gages  increase 
dramatically  in  both  tension  and  compression.  The  idealization  of  Fig. 
28  is  obtained  experimentally  only  if  the  strain  gages  are  located  at 
a  buckle.  This  implies  that  a  prior  knowledge  of  the  buckling  pattern 
is  necessary.  The  strain  gage  placement  for  all  30  panels  was  based 
upon  a  symmetrical  pattern  about  the  panel's  center  as  was  indicated 
by  STAGSC-1.  As  will  be  explained  later,  the  actual  pattern  did  not 
have  this  symmetry  which  meant  that  in  most  cases,  the  strain  gages  were 
not  located  at  a  buckle.  However,  good  results  were  still  obtained  even 
though  the  plots  of  strain  vs.  load  were  not  as  well  defined  as  in  Fig. 


The  least  reliable  method  of  determining  the  buckling  load  was  the 
Southwell  method.  The  Southwell  method  (see  Appendix  A  for  a  fundament¬ 
al  discussion)  required  that  the  horizontal  LVDT's  be  placed  in  very 
close  proximity  to  the  largest  buckle.  Since  a  matrix  made  up  of  nine 
and  at  times  twelve  LVDT's  was  used,  this  requirement  was  usually  met. 

On  a  number  of  panels  (see  Table  I)  this  method  gave  obviously  erroni- 
ous  results.  This  was  due  to  a  linear  relationship  of  applied  load  to 
radial  displacement.  Only  with  those  data  plots  that  exhibit  a  nonlin¬ 
ear  curve  could  a  reasonable  buckling  load  be  obtained.  The  linearity 
may  be  explained  by  the  combination  of  panel  imperfections  and  panel 
support  system  eccentricities  which  serve  to  strengthen  the  panel.  This 
prevented  the  relatively  larger  displacements  at  high  load  levels. 

In  order  to  measure  radial  eccentricities,  a  crude  measuring  device 
was  made  which  was  used  in  conjunction  with  a  set  of  dial  vernier  cal¬ 
ipers,  see  Fig.  29.  The  device  was  used  to  determine  panel  "out-of- 
roundness"  brought  about  by  less  than  perfect  simple  support  installa¬ 
tion.  Radial  measurements  were  taken  at  discrete  locations  on  11  panels. 
It  was  found  that  the  final  buckled  configuration  of  these  panels  gen¬ 
erally  correspond  to  the  pattern  of  their  respective  prebuckled  radial 
eccentricities. 

Buckling  Pattern 


It  was  found  that  the  installation  sequence  of  a  panel  into  the  test 


set-up  was  crucial  to  the  formation  of  its  buckling  pattern.  Fig's.  30 
to  41  show  the  buckling  pattern  progression  of  various  panels  as  com¬ 
pared  to  the  pattern  generated  by  STA6SC-1.  Notice  that  the  experi¬ 
mental  pattern  did  not  exhibit  the  symmetry  that  theory  would  suggest. 
In  fact,  when  the  vertical  supports  were  positioned  as  in  Fig.  12  with 
the  panel  protruding  i  of  an  inch  at  the  top  of  the  supports,  none  of 
the  panels  showed  any  buckling  in  their  lower  half.  This  is  probably 
due  to  the  panel's  inherent  flexibility  in  the  upper  half  brought  about 
by  the  unsupported  i  inch  space  on  each  side.  The  result  is  an  uneven 
stress  distribution  from  the  flexible  upper  half  to  the  stiffer  bottom 
half,  thus  preventing  buckling  at  the  bottom.  The  top  picture  in  Fig's 
31,  35,  36,  38,  and  40  shows  the  initial,  small  buckle  which  formed  at 
the  upper  right  portion  of  most  of  the  panels  tested.  This  buckle  man¬ 
ifested  itself  as  only  a  1  to  2  percent  drop-off  in  load,  or  a  kink  in 
the  load  vs  end-shortening  data.  As  the  load  was  increased,  the  buckle 
progressed  from  right  to  left  until  a  very  audible  and  visual  snap 
occurred  signifying  panel  collapse  and  a  distinct  drop-off  in  load. 

The  resulting  buckling  pattern  was  generally  observed  to  contain  three 
half-sine  waves  in  the  circumferential  direction  for  the  16  and  12  inch 
chord  length  panels  and  one  half-sine  wave  for  the  8  inch  chord  panels. 
This  is  compared  to  the  analytical  patterns  (Fig's.  32,  34,  37,  39,  and 
41)  which,  depending  upon  aspect  ratio  and  ply  orientation,  contain  3 
to  6  half-sine  waves  in  the  circumferential  direction  and  4  to  12  half¬ 
sine  waves  in  the  vertical  direction.  The  formation  of  the  initial 


Figure  30 


buckle  can  be  explained  by  eccentricities  in  the  test  set-up  since 
theory  indicates  that  the  first  buckle  should  occur  at  the  buckling 
load.  In  an  attempt  to  explain  this  deviance  from  theory,  the  panel 
installation  process  was  reexamined.  Initially,  all  the  panels  were 
installed  by  entirely  removing  the  left  simple  support.  The  right 
support  was  left  in  place  with  the  set  screws  untightened  to  separate 
the  knife  edges.  The  panel  was  then  slid  in  and  clamped  along  the  top 
and  bottom  curved  edges.  At  this  point  the  left  support  was  reinstalled 
and  set  screws  tightened  on  both  left  and  right  supports.  The  load  was 
applied  and  the  resulting  pattern  photographed. 

A  new  installation  procedure  was  followed  and  panels  1931,  1932, 
and  4931  were  retested.  This  time  both  left  and  right  supports  were 
removed.  This  enabled  the  panels  to  be  clamped  along  their  top  and 
bottom  edges  without  interference,  thus  preventing  a  possible  source  of 
panel  binding.  The  simple  supports  were  then  installed,  knife  edges 
meticulously  centered,  and  the  set  screws  finger  tightened.  This  pro¬ 
cedure  was  modified  somewhat  for  the  1932  panel.  The  i  inch  free 
space  was  placed  at  the  bottom  corners  of  this  panel.  This  was  accom¬ 
plished  by  lifting  the  knife  edges  until  they  contacted  the  upper  plate 
located  at  the  top  horizontal  cross-head.  The  set  screws  were  tightened 
just  enough  to  keep  the  knife  edges  from  sliding  back  down. 

The  new  installation  procedures  produced  significant  results.  It 
was  observed  that  the  final  buckled  configuration  of  both  the  1931  and 


4931  panels  was  the  same  as  before  the  new  Installation  procedures  were 
implemented,  see  Fig's.  42  and  43.  However,  there  was  no  buckling  pro¬ 
gression,  indicating  that  the  new  installation  procedure  corrected  the 
previous  eccentricity  problems.  As  per  theory,  both  panels  maintained 
their  original  shapes  until  the  critical  buckling  load  was  reached;  at 
which  time  they  instantaneously  snapped  into  their  buckled  configura¬ 
tions.  It  is  important  to  note  that  the  respective  buckling  loads 
associated  with  the  two  different  installation  procedures  differed  by 
only  3%.  The  most  dramatic  results  were  obtained  when  the  i  inch  free 
panel  space  was  placed  at  the  bottom.  This  configuration  enabled  the 
1932  panel  to  buckle  into  the  first  symmetrical  pattern  observed  in 
this  thesis,  see  Fig.  44.  It  would  appear  that  when  simple  supports 
are  used,  the  panel's  free  space  must  be  on  the  same  end  as  the  applied 
load.  This  allows  for  a  more  uniform  stress  distribution  which  results 
in  buckling  symmetry. 

Ply  Orientation/Aspect  Ratio  Effects 

Given  a  specific  ply  orientation.  Table  II  shows  the  analytical 
bifurcation  loads  decreased  with  increasing  aspect  ratio.  This 
occurred  experimentally  with  the  (90,±45,0)g  and  (±45)2g  panels.  The 
(90,0)25  panels  showed  one  deviation  from  this  order  with  the  1-aspect 
ratio  panels  having  the  highest  buckling  load.  Analytical  values  also 


show  that  for  a  specific  aspect  ratio,  the  highest  buckling  load  was 
associated  with  the  (90,145,0)^  ply  orientation  followed  by  the  (±45)25 
and  (90,0)25  orientations,  respectively.  Experimentally,  this  sequence 
was  followed  by  only  the  i  and  3/4  aspect  ratio  panels.  (See  Appendix 
C  for  the  extentional  and  bending  stiffness  matrices  for  each  ply 
orientation. ) 

It  was  determined  that  the  aspect  ratio  had  the  greatest  effect  on 
a  panel's  buckling  load,  see  Fig.  45.  For  the  same  ply  orientation, 
the  maximum  Increase  In  aspect  ratio  (i  to  2)  produced  over  a  50%  drop 
in  buckling  load  both  analytically  and  experimentally. 
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Experimental  Buckling  Load 
vs.  Aspect  Ratio 

VI.  CONCLUSIONS 


Based  on  the  analytical  and  experimental  results  the  following 
conclusions  can  be  made  for  composite  cylindrical  panels  subjected 
to  axial  compressive  loading. 

1.  Though  the  author  sees  imperfections  as  a  very  important 
contributor  to  bifurcation  load,  this  thesis  deals  only 
with  the  linear  bifurcation  included  within  STAGSC-1. 

It  can  be  seen  that  a  linear  analysis  gives  good  results 
for  the  support  device  used. 

2.  The  variance  in  knockdown  factors  precluded  any  general 
correlation  between  it  and  a  panel's  ply  orientation  and/ 
or  aspect  ratio. 

3.  In  a  great  majority  of  cases  the  observed  buckling  pattern 
did  not  exhibit  the  symmetry  that  would  be  expected  by 
purely  analytical  means.  It  was  found  that  the  existing 
simple  support  device  was  the  cause.  Due  to  stress  dis¬ 
tribution  considerations,  the  i  inch  free  panel  space 
must  be  on  the  same  end  as  the  applied  load  for  a  sym¬ 
metrical  buckling  configuration  to  occur. 

4.  With  any  test  device  dealing  with  curved  panels,  the 
researcher  should  be  particularly  careful  of  the  panel 
boundary  restraints.  A  less  than  ideal  installation 
procedure  will  cause  panel  eccentricities  resulting  in 


initial  buckles  forming  which  do  not  correspond  to  the 
final  buckling  load  of  the  panel. 

5.  The  buckling  pattern  was  more  sensitive  to  boundary 
restraint  irregularities  than  was  the  buckling  load 
which  remained  essentially  unchanged  for  a  given  panel. 

6.  It  appears  that  the  STA6SC-1  significantly  overestimates 
the  number  of  half-sine  waves  in  the  panels  buckled 
configuration  when  compared  to  the  number  occurring 
experimentally. 

7.  The  strain  reversal  and  Southwell  methods  for  determin¬ 
ing  buckling  are  generally  inferior  to  the  end-shortening 
method  because  strain  reversal  requires  a  prior  knowledge 
of  the  buckling  pattern  and  the  Southwell  method  is  par¬ 
ticularly  sensitive  to  panel  eccentricities.  However, 
from  Table  II,  it  can  be  seen  that  they  provide  a  good 
check  of  the  end-shortening  method. 

8.  The  amount  of  load  a  panel  can  carry  is  dependent  upon 
its  aspect  ratio  and  ply  orientation.  It  was  found  that 
for  the  boundary  conditions  used  in  this  thesis,  the 
aspect  ratio  had  the  greatest  effect  upon  the  buckling 
load  followed  by  the  ply  orientation. 

9.  During  the  actual  experimental  loading  of  the  panels, 

it  was  observed  that  the  (90, ±45,0).  and  (90,0),.  panels 


were  more  prone  to  breakage  and  splintering  at  the 
top  corners  than  were  the  (±45)2^  panels.  This  appears 
to  be  due  to  the  overall  diagonal  stiffening  character¬ 
istic  of  the  45®  oriented  fibers. 

10.  As  can  be  seen  in  Fig.  45,  the  difference  in  buckling 
loads  associated  with  the  three  ply  orientations  used 
in  this  thesis  generally  decrease  with  increasing  aspect 
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The  Southwell  Method 


The  Southwell  Method  [19]  is  a  technique  for  obtaining  the 
buckling  load  of  a  structure  from  experimental  load  deflection  infor¬ 
mation.  It  has  been  used  for  the  buckling  of  columns,  beam  columns, 
plates,  and  more  recently,  shells  [5].  The  method  involves  using  the 
well  established  (free  of  test  fixture  slack)  portion  of  an  applied 
load  (P)  vs  lateral  deflection  (fi)  curve.  Selected  points  are  then 
normalized  to  obtain  a  plot  of  6/p  vs  6.  A  linear  least  squares  fit 
of  the  points  gives  the  slope  of  the  line 

6/P  _  1 

6  ■  P  (1 

Inverting  the  slope  gives  the  experimental  buckling  load,  P.  Figure 
1-A  shows  P  vs  6  curve  for  the  3911  panel .  The  curve  appears  to  be 
well  established  between  3000  and  5000  lbs.  Data  reduction  was  per¬ 
formed  on  9  points  in  the  above  load  range  and  is  shown  on  Fig.  2-A. 
The  Southwell  buckling  load  was  determined  to  be  5996.7  lbs. 
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Figure  1-A  Load  vs.  Radial  Displacement  for  3911  Panel 


Appendix  B 


STAGSC-1  Examples 

T1m6/$pace  Limits  The  amount  of  computer  time  and  space  required 
by  STAGSC-1  is  dependent  upon  a  panel's  number  of  degrees  of  freedom. 

In  this  thesis,  the  16x16  inch  panels  contained  the  most  active  degrees 
of  freedom  with  5985,  while  the  8x16  inch  panels  contained  the  least, 
2913.  The  following  maximum  limits  were  found  to  be  sufficient  for  all 
cases:  CP  Time  «  1200  seconds,  10  Time  -  5000  seconds,  and  CM  -  176,001 
words.  Figure  1-B  shows  the  input  data  deck  for  the  (±45)2^,  16"xl6" 
panel.  Figure  2-B  shows  the  necessary  control  cards  for  the  program. 
Consult  reference  [14]  for  further  information. 
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Figure  1-B  Input  Data  Deck  for  the  Linear 

Bifurcation  Analysis  of  the  193  Panels 


1.  HEB,  T1200,  I/O  5000,  CM  176000.  D710554. 

2.  AHACH,  PROCFIL,  ID=A810171.  SN=ASDAD. 

3.  BEGIN,  NOSFILE 

4.  ATTACH.  STAGS!,  1D«D820138.  MR-1. 

5.  GET.  SIM45A.  ID-HEBERT. 

6.  STAGSl.  SIM45A. 

7.  RETURN,  STAGSl. 

8.  AHACH,  STAGS2.  ID-D820138. 

9.  STAGS2,  MR-1. 


Figure  2-B  Control  Cards  for  the  Linear  Bifurcation  Analysl 

of  the  193  Panels 
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Appendix  C 

Extenslonal  and  Bending  Stiffness  Matrices 


Extensional  Stiffnesses 


A 

A 

A 

11 

12 

16 

A 

A 

A 

12 

22 

26 

A 

A 

A 

16 

26 

66 

ti 


(90.±45,0)25 


.35 

.11 

0 


.11 

.35 
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0 

0 

.12 


(±45)25 

.26 

.20 

0 

.20 

.26 

0 

X 

0 

0 

.21 

X  10' 


(90.0)25 


10®  Ibf/ln 


Ibf/ln 


89 


Bending  Stiffnesses  (Cont'd.) 


39.12  2.34  0 

2.34  77.72  0 

0  0  4.00 


Ibf-ln 


Note:  All  coupling  tenns  (B^j)  are  zero 
since  the  three  ply  orientations 
are  synmetrlc  about  their  mid-surface. 

As  was  previously  mentioned,  the  buckling  of  shells  Is  basically  a 
bending  problem.  This  obviously  Implies  that  the  bending  stiffness 
matrix  Is  of  primary  Importance  In  determining  a  shell's  buckling  load. 
Since  an  explicit  equation  for  solving  for  the  buckling  load  of  a 
cylindrical  shell  was  not  available  to  the  author,  Eq.  (1C),  a  buckling 
equation  for  a  flat  composite  plate  with  an  aspect  ratio  of  2  [10],  was 
used  to  compare  the  relative  sizes  of  the  buckling  loads  associated 
with  the  3  ply  orientations  used  in  this  thesis. 

D 

N,  «  +  2(D  +  2D  )  +  4D  (1C) 

X  tt  12  66  22 

where  represents  a  force  per  unit  length.  Upon  substituting  the 
appropriate  values  for  the  Into  Eq.  (IC),  It  was  found  that  the 
(90,±45,0)9.  had  the  largest  buckling  load  followed  by  the  (90,0)9^ 


and  (±45)^,  respectively.  This  was  not  what  occurred  analytically  nor 
in  most  cases  experimentally,  where  it  was  found  that  the  (90,0)2^ 
orientation  buckled  at  the  lowest  load.  It  is,  therefore,  evident 
that  in  shell  analysis  the  extensional  stiffness  terms  (A.^)  must 
also  be  included. 
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